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Abstract 

We show that both the muhiphcity distribution and the ratio of factorial 
cumulants over factorial moments for 2-jet events in e+e~ annihilation at the 
Z° peak can be well reproduced by the weighted superposition of two negative 
binomial distributions, associated to the contribution of hb and light flavoured 
events respectively. The negative binomial distribution is then suggested to 
describe the multiphcity distribution of 2-jet events with fixed fiavour. 
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1. Introduction 



Two different experimental effects in the Multiplicity Distributions (MD's) of charged 
particles in full phase space in e^e~ annihilation at the Z° peak, i.e., the shoulder visible 
in the intermediate multiplicity range [0, ^, |^ and the quasi-oscillatory behaviour of the 
ratio of factorial cumulants over factorial moments of the MD, Hq, when plotted as a 
function of its order g ^, have been quantitatively reproduced in |]^ in terms of a 
weighted superposition of two Negative Binomial Distributions (NBD's), associated to 
two- and multi-jet events respectively. A further test of this picture, in which the simple 
NED appears at a very elementary level of investigation, is provided by the study of 
samples of events with a fixed number of jets. In [0, the Delphi Collaboration has 
shown that a single NBD can describe the MD's for events with a fixed number of jets for 
a range of values of the jet resolution parameter j/min- This was indeed the starting point 
of the parametrization proposed in ||^. 

In this letter, by extracting the ratio Hq from published MD's according to the pro- 
cedure described in [^, we show that the oscillations observed experimentally are larger 
than those predicted by a single NBD, even after taking into account the truncation ef- 
fect, which was shown Q to be important in the behavior of Hq's. These results suggest 
that, while hard gluon radiation plays a relevant role in the explanation of the shoulder 
structure of MD's and of oscillations of the ratio Hq, some other effects should be taken 
into account for a detailed description of experimental data of events with a fixed number 
of jets. In this respect, it is worth recalling the interesting results obtained by the OPAL 
Collaboration on forward-backward correlations and on the increase of transverse mo- 
mentum of produced hadrons in the intermediate multiplicity range: it has been found 
indeed that both effects are mainly due to hard gluon radiation, i.e., to the superposition 
of events with a fixed number of jets. However, a residual positive correlation has been 
found in a sample of 2-jet events; via Monte Carlo simulations, this effect has been as- 
sociated to the combined action of superposition of events with different quark flavours 
and, in the central region, to a residual effect due to resonances' decays. Let us remind, 
however, that the presence of heavy flavours has been shown not to affect the increase 
of the transverse momentum of produced hadrons in the intermediate multiplicity range, 
thus suggesting that not all observables are sensitive to the original quark flavour. A 
theoretical study based on Monte Carlo simulations has first suggested that the study of 
MD's can indeed point out interesting features of particle production in hb event s |lTO[] . Re- 



cently the Delphi Collaboration has established experimentally the sensitivity of MD's 
to the original quark flavour, by comparing the MD for the full sample of events with a 
sample enriched in hh events [0]. 

We propose in this letter to associate a NBD to the MD in 2-jet events of fixed flavour. 
We show, after examining possible alternatives, that the weighted superposition of two 
NBD's, which we associate to hh and light flavoured events, describes very well both the 
MD's and the ratio Hq for 2-jet events; the two NBD's have the same k parameters and 
differ in the average multiplicity only. Some consequences of this fact are examined in the 
conclusions. 
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Table 1: Parameters and x^/NDF of the fit to experimental data[|ril] on single hemisphere 
MD's for bb sample and for all flavours with a single NBD and with a NBD shifted by one 
or two units. 





bb sample 


all flavours 


NBD 


n 


11.67±0.07 


10.67±0.02 


k 


24±2 


14.5±0.3 


xVndf 


118/26 


140/28 


NBD (shi 


"t by 1 unit) 


n 


10.62±0.07 


9.63±0.02 


k 


15.8±0.6 


10.42±0.2 


xVndf 


60/26 


64/28 


NBD (shi 


"t by 2 units) 


n 


9.55±0.07 




k 


10.4±0.3 




xVndf 


20/26 





2. MD's in fe-jets 

The Delphi Collaboration has studied the effect of quark flavour composition on the MD 
in one hemisphere, by comparing the MD for the full sample of events with that for a 
sample enriched in bb events [pi] ]: the MD extracted from the bb sample was found to be 
essentially identical in shape to the MD obtained for the full sample, apart from a shift 
of one unit, which may be related to the effect of weak decays of i?-hadrons[0. To give 



a quantitative comparison of the MD's in single hemisphere in the bb sample and in the 
sample with a mixture of all flavours, we have fitted both experimental MD's with a NBD 
and with a NBD shifted by one or two units. The results of the fit are shown in Table |]. A 
single NBD gives a poor description of the MD for the bb sample; the description improves 
strongly if one introduces a shift by one unit, and becomes even better with a shift by 
two units. The reason is that with the shift the NBD is able to reproduce better the 
head of the distribution; however the tail remains underestimated. In this respect, one 
should remember that the single NBD cannot give a good description of the MD with all 
flavours in full phase space, since it cannot reproduce the shoulder structure due to the 
superposition of events with different number of jets|0]. The fact that this feature should 
be present for the hb sample too should be verified experimentally. 

In any case, interesting information on the properties of these MD's can be extracted 
without using any parametrization at all. In what follows we will consider only 2-jet 
events, selected with a suitable algorithm, but the same reasoning can be carried out also 
for the full sample. Let us call p„, and the MD's in a single hemisphere for all events, 
for bb events and for light flavoured (non bb) events respectively, and g{z) = J2'^=oPnZ"', 
g^{z) and g^z) the associated generating functions. With a the fraction of bb events, one 
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has: 



I.e., 



Pn = oip^ + (1 - 



g{z) = ag\z) + (1 - a)g\z) (2) 

These relations are vahd in general; Delphi data and our analysis shown in Table |l] 
suggest that is given by p„ with a shift of one unit: 



Pn = Pn-l 



n > 0; 



P'o 



0; 



I.e., 



Substituting now eq. (11) in eq. 



g'^iz) = zg{z) 
]), one gets: 

I — a , , 



I.e., 



9\z) 



az 



-9[z) 



rz(l — a) 



az 



9\z) 



(3) 



(4) 



(5) 



(6) 



The MD in hh events is the convolution of a shifted geometric distribution, of average 
value 1/(1 — a), with the MD in light flavoured events. The shifted geometric MD could 



be related to the MD of the decay products of B hadrons in the framework of |T2|. 

The connection of the MD in a single hemisphere to that in full phase space is not 
entirely trivial, since one has to take into account additional effects, like for instance 
charge conservation, which requires that the final multiplicity be even. 

Let us use the same notation for MD's as in the previous paragraph, but with capital 
letters to denote the MD's in full phase space; by taking the two hemispheres as inde- 
pendent (which, as suggested in 0, should be a good approximation at least for light 
flavours) but applying charge conservation, we obtain: 



P\n^,n2) 







if Til + n2 is even 
otherwise 



(7) 



Here P^{ni, is the probability to produce rii particles in one hemisphere and n2 in the 
other hemisphere and i denotes either all 2-jet events (no label), or hb events (z = h) or 
light flavoured events [i = I). The factor 2 is for normalization, assuming that the do 
not privilege the even or odd component; in any case this effect can be easily taken into 
account and results similar to those given below are obtained. 
The MD in full phase space is given by definition by: 



n 



ni) 



ni=0 



In terms of the generating functions, one obtains 



00 n 



G\z) = ^""Pn = J2 J2 z'''p\ni)z''-"'p\n - n^) 



n=0 



n=0 ni=0 
n even 



m^)?+[9\-z)r) 



(9) 



4 



We can see now how the relations which we obtained in the previous paragraph are 
modified going from a single hemisphere to full phase space: by putting eq. (H) in eq. (H), 
one has 

For a = 0, one obtains that in full phase space the MD in bb events coincides with the 
MD for light fiavoured events with a shift of two units, as we get G*(-2) = z'^G\z). The 
MD for small values of a should not be too far from this limit, as can be easily checked 
with numerical examples. 

In conclusion, going from single hemispheres to full phase space by taking into account 
only charge conservation, the MD of bb events becomes close not to the total MD but to 
the MD of light flavoured events. The two MD's seem to have the same characteristics, 
like average value and dispersion, and the only difference should lie in a shift of two units. 



3. MD's and HqS ratio in 2-jet events 

The analysis of MD's for events with a fixed number of jets, and in particular 2-jet events, 
has been performed in [0, where a single NBD has been shown to reasonably describe 
the MD's for events with a fixed number of jets, for several values of the jet resolution 
parameter ?/niin (the JADE jet-finding algorithm has been used in 0). A comparison of 
Delphi data with yuan = 0.02 with a single NBD is shown in Figure |l|a together with the 
residuals, i.e., the normalized difference between data and theoretical predictions, which 
point out the presence of substructures in experimental data. In view of previous results, 
it is then interesting to investigate whether these substructures can be explained in terms 
of the different contribution of quark-jets with different flavours. 

We parametrize the experimental data on MD's for 2-jet events in full phase space in 
e^e~ annihilation at the peak[0 in terms of the superposition of 2 NBD's, associated 
to the contribution of b- and light flavours (we include the charm among the light flavours 
to a flrst extent). We flx therefore the weight parameter to be equal to the fraction of 
bb events, a = 0.22||T^. Following the results of the previous section, we ask that the 
NBD associated to the b flavour be shifted by two units and that both parameters of the 
two NBD's, n and k, be the same. Formally, we perform then a flt with the following 
2-parameter distribution: 

Pnin, k) = aP^^^in, k) + {1 - a)P^''{n, k) (11) 

where P^2(^5 fc) = for n < 2. P^^{n, k) is here the NBD, expressed in terms of two 
parameters, the average multiplicity n and the parameter k, linked to the dispersion by 
£)2/^2 = i/n^i/k^ as: 

,NB/- 7\ k{k + 1) . . .{k + n -1) ( k \^ f n 



As far as MD's in full phase space are concerned, one has also to take care of the "even- 
odd" effect, i.e., of the fact that the total number of flnal charged particles must be even 
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Figure 1: a): charged particles' MD for two-jet events in full phase space, P„, at the 
peak from Delphi[|^ with y-aan = 0.02 are compared with the fit with a single NBD 
as performed by Delphi Collaboration (solid lines); the lower part of the figure shows 
the residuals, i.e., the difference between data and theoretical predictions, expressed 
in units of standard deviations. The even-odd effect has been taken into account (see 
eq. (0)). b): Same as in a), but the solid line here shows the result of fitting eq. (|llD, 
with parameters given in Table |^. c) : Same as in a) , but the solid line here shows the 
result of fitting eq. ([l^) , with parameters given in Table ^. 
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Table 2: Parameters and x^/NDF of the fit to experimental data on 2-jet events MD's 
from Delphi [0 with three different MD's: the weighted superposition of a NBD and a 
shifted NBD with the same parameters (eq. ([TTD ), the weighted superposition of a Poisson 
plus a shifted Poisson (eq. (0)) and the weighted superposition of two NBD's with the 
same k (eq. (|20|) ). The weight used is the fraction of bb events. The even-odd effect 
has been taken into account (see eq. (|13|)). Results are shown for different values of the 
jet-finder parameter j/min- 





Vmin = 0.01 


ymin = 0.02 


Vmin = 0.04 


NBD + shifted NBD (same n and A;) eq. (11) 


n 


17.17±0.05 


18.01±0.04 


18.99±0.05 


k 


69±5 


57±3 


44±1 


XVNDF 


18.2/17 


27.9/17 


53.9/21 


Poisson + shifted Poisson eq. (|14|) 


rii 


19.67±0.13 


21.10±0.10 


22.83±0.09 


n2 


16.37±0.07 


16.95±0.06 


17.51±0.06 


XVNDF 


21.7/17 


20.75/17 


45.05/21 


2 NBD (same k) eq. (0) 


ni 


16.81±0.21 


17.22±0.15 


17.98±0.15 


fib 


20.26±1.71 


21.96±1.57 


23.61±1.64 


k 


124±51 


145±53 


120±33 


XVNDF 


17.4/16 
3.44±0.83 


12.6/16 
4.6±0.5 


27.5/20 
5.6±0.5 



due to charge conservation; accordingly, the actual form used in the fit procedure is given 
by: 

pfps ^ I APn if n is even ^-^g-^ 
" 1 otherwise 

where A is the normalization parameter, so that Zl^o-^rf^* ~ ^■ 

Table |^ (first group) shows the two parameters n and k of eq. ( pi]) for different values 
of the jet resolution parameter j/min; the proposed parametrization gives a rather good 
description with only two parameters; the agreement is worse for ymin = 0.04, which 
could be due to the contamination of 3-jet events. However, as shown in Figure |l|b 
for ?/min = 0.02, the oscillatory structure in the residuals does not disappear with this 
parametrization. Let us remind that the values of x^/NDF shown in the Table should 
be considered just indicative, since we did not know the full covariance matrix and we 
could not then treat properly the correlations between different channels of the MD. This 
forbids also a direct comparison of the x^/NDF of the present parametrization with the 
values of x^/NDF for a single NBD obtained in where correlations between bins were 
taken into account. Finally, let us say that we also fit eq. (P) by assuming that pl^ is a 
NBD, and we found the same results. 

It is interesting at this point to investigate a minimal model where no physical corre- 
lations are present at the level of events with fixed number of jets and fixed flavour: we 
have performed a fit using a weighted sum of a shifted Poisson plus a Poisson distribution 
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(with the correction for the even-odd effect according to eq. (f3|)): 



Pn{ni,n,) = aPl_^{n,) + (1 - a)pP(n2) (14) 

where 

P^{n) = -e-^ (15) 
nl 

Also in this case, we have two free parameters. Resuhs of the fit are shown in Table |] 
(second group); also in this case a reasonable fit is achieved, even though the MD at 
1/min = 0.04 shows again some anomaly. The MD with the parameters shown in Table |] 
is compared to experimental data with = 0.02 in Figure [^c; it should be pointed out 
that the structure in the residuals is present in this case too. From this analysis, one 
would then conclude that physical correlations visible in e+e" annihilation result trivially 
from the superposition of samples of events with different quark flavours. A more accurate 
analysis with full covariance matrix is needed to see which parametrization is preferred 
by experimental data. However, independent of the chosen parametrization, two different 
components, which can be associated to b- and light flavours contributions, are visible in 
the MD of 2-jet events. 

A more detailed analysis of the tail of the MD, which can help distinguish different 
parametrizations, comes from the study of the ratio of unnnormalized factorial cumulant 
over factorial moments 

H, = ^ (16) 

as a function of the order q. The factorial moments, Fq, and factorial cumulant moments, 
Kg, can be obtained from the MD, P„, through the relations: 

oo 

Fg = Y.n{n-l)...{n-q + 1)P„, (17) 

n=q 

and ^ 

^, = F,-Y.[^'^^^Kq_,h (18) 

Since the Hq's were shown to be sensitive to the truncation of the tail due to the finite 
statistics of data samples^, moments have to be actually extracted from a truncated MD 
defined as follows (including again the correction for the even-odd effect as in eq. (p!3D): 



ytrunc 



A'Pn if {rimin <n< Umax) and n is even ^^g^ 
otherwise 



Here n^m and Umax are the minimum and the maximum observed multiplicity and A' is 
a new normalization parameter. 

In Figure ^ the H^s extracted from the experimental MD published by Delphi Col- 
laboration (here ymin = 0.02) with the procedure explained in are compared with 
the predictions of a single NBD as fitted by Delphi Collaboration []^, and of eq.s (|TT|) 
and (|T^). It is clear that all three parametrizations fail to describe the experimental be- 
haviour of the ratios Hq, i.e., the description of the tail of the MD is not accurate. We 
then conclude that a single NBD cannot describe accurately the MD in 2-jet events, as 
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already suggested by the study of residuals of MD's. The superposition of two NBD's 
with the same parameters turns out also to be inadequate; one concludes that the im- 
posed constraints are too strong and that some additional differences between bb and light 
flavoured events should be allowed. Finally, the observed deviation from the Poisson-like 
fit suggests that there are indeed dynamical correlations beyond the purely statistical 
ones. 



4. A new parametrization of MD's in 2-jet events 

In the previous paragraph, the difference between the average multiplicity in bb and light 
flavoured events in the parametrization (pUf ) has been fixed to 2; by using the parametriza- 
tion dl^), where this difference is not a priori constrained, larger values have been obtained. 
Let us also remind that the experimental value of this observable at the peak is close to 
2.8|T^; theoretical predictions in the framework of Modified Leading Log Approximation 



plus Local Parton Hadron Duality |]T5[ give even larger values. It is therefore interesting to 
investigate whether one can reproduce not only the shape of the MD, but also its tail and 
then the ratio Hg, by using a superposition of two NBD's, but relaxing the constraint on 
the average multiplicities. The only constraint we impose is that the parameters k of the 
two NBD's be the same, while we allow a variation of the difference between the average 
multiplicities. For the sake of simplicity and in order to be more independent from any 
theoretical prejudice, we do not include any shift in the MD for bb events. 

Formally, we perform then a fit with the following 3-parameter MD (plus the correction 
for the even-odd effect in eq. (P^)): 

P^ni, Uh, k) = aP^^fi,, k) + {l- a)P^''{nu k) (20) 

The parameters of the fits and the corresponding x^/NDF are given in Table ^ (third 
group) for different values of the resolution parameter 2/min- A really accurate description 
of experimental data is achieved. Notice that the best-fit value for the difference between 
the average multiplicities in the two samples, 5hu also given in Table |^, is quite large. This 
difference grows with increasing ?/min, i-e., with increasing contamination of 3-jet events. 

Figure (|^) compares the predictions of eq. ( PP] ) with the experimental MD's for two-jet 
events at different values of the resolution parameter yuan- The residuals are also shown 
in units of standard deviations. One concludes that the proposed parametrization can 
reproduce the experimental data on MD's very well; no structure is visible in the residuals. 

As already discussed, the ratio Hg gives a more stringent test of theoretical parametriza- 
tions; it is then interesting to study the predictions of eq. ( pOD for this ratio. In this case, 
one can obtain a closed expression for the factorial moments in terms of the parameters 
5bu ni and k. Let us notice indeed that, since the two components are given by a NBD 
with the same /c, they have the same normalized factorial moments, which for a NBD are 
given by: 

^(0 = Fj^) = + ^) (21) 
i=i ^ 

From eq. (pOf), one obtains a similar relation for the generating function: 

G{z) = aG^''\z) + {I - a)G'^'\z) (22) 
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Figure 2: The ratio of factorial cumulant over factorial moments, Hg, as a function of q. 
a): Experimental data (diamonds) for 2-jet events with ?/min = 0.02 are compared with 
the fit with a single NBD as performed by Delphi Collaboration (solid lines), b): Same 
as in a) , but the solid line here shows the result of fitting eq. (pUD , with parameters given 
in Table 0. c): Same as in a), but the solid line here shows the result of fitting eq. (p!^, 
with parameters given in Table ^ In all three cases the even-odd and the truncation 
effects have been taken into account (see eq. (P!^)). 
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Figure 3: Charged particles' MD for two-jet events in full phase space, P„, at the 
peak from DELPHI []^ with different values of y^am are compared with a fit with the sum 
of 2 NBD's with the same parameter k as in eq. (^Ol). The even-odd effect has been taken 
into account (see eq. ([l3|) ). The lower part of the figure shows the residuals, Rn, i.e., 
the difference between data and theoretical predictions, expressed in units of standard 
deviations. 
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By differentiating the previous equation, one then gets the following expression for the 



unnormalized factorial moments, Fg-. 



n 



Fi=ni + a5t 



(23) 



(0 



aFf + (1 - 



oi{ni + + (1 - a)nl 



(24) 



9-1 



i=l 



n 1+ 



k 



Predictions of the ratio Hq as a function of the order q are obtained by inserting 
eq. (|20| ) into eq. (|19|). These predictions with parameters fitted to reproduce the MD as 
given in Table |^ (third group) are compared in Figure § with the HqS extracted from 
experimental data on MD's for 2-jet events0 at different values of ?/min according to the 
procedure described in [^]. The new parametrization gives an accurate description of the 
shape of MD's and is shown to describe well the ratio Hq too. Small deviations are still 
present for ?/min = 0.01, where 2-jet events are more collimated. They might be due to 
more subtle not yet understood effects. This consideration notwithstanding, the overall 
description of 2-jet MD's and HqS appears quite satisfactory. This result gives therefore 
further support to the parametrization of the MD in quark-jets with fixed flavour in terms 
of a single NBD. It is also remarkable that the average number of particles only depends 
on flavour quantum numbers, whereas the NBD parameter k is flavour-independent. 

As a further check, we also investigated the MD of 2-jet events with fixed flavour in the 
Monte Carlo program Jetset 7.4 PS|jl6[. For each flavour, we generated a sample of 60000 
events and of 60000 2-jet events (selected using the JADE algorithm with ?/min = 0.02) by 



using the OPAL tuning[jTj] and we fitted the MD's in full phase space with a single NBD, 
eventually with a finite shift. In the all-events sample, the x^/NDF is really bad, thus 
indicating that a single NBD cannot describe the MD of events with fixed flavour. By 
requiring the 2-jet selection, the description improves strongly; the MD for light quarks 
are indeed well reproduced by a single NBD, while hb events are better described by a 
shifted NBD. 



5. Conclusions 

It has been shown that a single NBD cannot reproduce the observed behavior of the ratio 
Hq for events with a fixed number of jets in full phase in e^e~ annihilation at the 
peak. A simple phenomenological parametrization of the MD in terms of the weighted 
superposition of two NBD's has been shown to describe simultaneously both the MD's 
and the ratio Hq. The weight of the first component was taken equal to the fraction 
of hh events, i.e., the two components were identified with the h- and the light flavours 
contribution respectively. The simple NBD parametrization is thus reestablished at the 
level of 2-jet events with fixed quark flavour composition. It is interesting to note that 
this result is consistent with the results obtained in the context of a thermodynamical 



model of hadronizationflS 



It is remarkable that the two NBD's associated to h- and light flavours contributions 
have the same parameter /c; since /c~^ is the second order normalized factorial cumulant, 
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Figure 4: The ratio of factorial cumulant over factorial moments, Hq, as a function of g; 
experimental data (diamonds) for 2-jet events with different values of y-cmn are compared 
with equation (^) (solid lines). The parameters used are shown in Table |^. The even-odd 
and the truncation effects have been taken into account (see eq. (|T^). 



13 



i.e., it is directly related to two-particle correlations, one concludes that two-particle 
correlations are flavour-independent in this approach. In addition, since both MD's are 
well described by a NBD, higher order correlations show a hierarchical structure [|1^], which 
is also flavour independent. This result can also be interpreted in the framework of clan 
structure analysis ||2CI||, where gets the meaning of an aggregation coefficient, being 
the ratio of the probability to have two particles in the same clan over the probability 
to have the two particles in two separate clans: in this language, one concludes that the 
aggregation among particles produced into clans in bb and light flavoured events turns 
out to be the same. The flavour quantum numbers affect then the average multiplicity 
of the corresponding jets only, but not the structure of particle correlations. It would 
be interesting to see, when appropriate samples of events will be available, whether this 
property established in full phase space continues to be valid in restricted regions of it. 
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